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Abstract. This paper discusses a recent proposal for the simulation of acoustic black 
holes with ions [1]. The ions are rotating on a ring with an inhomogeneous, but 
stationary velocity profile. Phonons cannot leave a region, in which the ion velocity 
exceeds the group velocity of the phonons, as light cannot escape from a black hole. 
The system is described by a discrete field theory with a nonlinear dispersion relation. 
Hawking radiation is emitted by this acoustic black hole, generating entanglement 
between the inside and the outside of the black hole. We study schemes to detect the 
Hawking effect in this setup. 
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Figure 1. Schematic depiction of the ion rotation with velocity v (9) and the 
phononic excitations with velocity ±c(#), that depend on the varying ion spacing. 
In the subsonic region the ion velocity is v (9) = v m i n , in the supersonic region it 
is v (9) — w max . The external force F° (9) localized at the transition between the 
super- and subsonic regions (de-) accelerates the ions. The black area represents the 
supersonic black hole region. A movie of the ion motion is accessible online (New J. 
Phys. 13, 045008 (2011)). 

1. Introduction 

It was shown based on the theory of quantum fields (QFT) in curved spacetime that, 
surprisingly, black holes emit thermal radiation [2]. However, the direct observation 
of Hawking radiation is difficult because the Hawking temperature is very small for 
astrophysical black holes. Furthermore, the original derivation of Hawking radiation 
relies on the validity of the wave equation on all scales, whereas QFT in curved space 
is expected to be reliable just up to the Planck scale. This paper is discussing Hawking 
radiation in the context of analogous hydrodynamical systems, which address these 
issues. These have an analog horizon at the transition from subsonic to supersonic (black 
hole) flow. It was shown that such analog horizons emit Hawking radiation [3], which 
can potentially be detected in experiments. The variety of different analog models and 
their tunability allows us to study the robustness of Hawking radiation against changes 
in the underlying microphysics. This will contribute to deepen our understanding of the 
Hawking effect also in gravitational black holes. 

In recent years many experimental tests of Hawking radiation have been proposed 
based on the hydrodynamical analogy [I], e.g., phonons in Bose-Einstein condensates 
(BECs) [51 [6], Fermi gases [7], superfluid Helium [8], slow light P [101 EH]) an d 
nonlinear electromagnetic waveguides [12J. Some experiments have implemented analog 
spacetimes in experiments, e.g., in BECs [13], optical fibres [II] . and water surfaces [T5]. 
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Figure 2. Typical profiles of the ion [v, green) and the sound velocity (c, blue) used 
in this article (see | Appendix B[ ). The profile shows a super- and a subsonic region. The 
ion velocity in the subsonic region is v m \ n = 2ir x 0.83/T and the black hole horizon is 
located at <Jv m - m T = 2ir x 0.25. A white hole horizon is also present on the ring pQ. 



Common to all these attempts is that they have not reported on the observation of the 
quantum Hawking effect [j]. 

In the present work we discuss a proposal how to build an analog model of a black 
hole in an experimentally realizable system of ions [lj. A special ingredient of our 
proposal is its discreteness, which naturally leads to a sublinear dispersion relation at 
high wavenumbers. This affects the trajectories of blue shifted waves close to the event 
horizon [IT]. The dispersion relation is, additionally, non-trivial at low wavenumbers 
because of the long range Coulomb force. However, as we will show, we still obtain 
Hawking radiation. Analytic derivations [TFJ HHJ HH [20] as well as numerical calculations 
show that the Hawking effect is robust against such short scale modifications, e.g., 
for a continuous field with a sublinear dispersion relation [21] and a discretized field 
on a falling lattice [IB]. Our proposal uses a parameter regime which is accessible in 
experiments at temperatures currently achieved. Thus, it could lead to the experimental 
observation of Hawking radiation. 

Let us now summarize the main idea of our proposal. We are constructing a discrete 
analog of an hydrodynamical system with super- and subsonic regions in a quadrupole 
ring trap [221 [23] as schematically depicted in Fig. [U The ions are rotating on a ring 
with circumference L with an inhomogeneous velocity profile v (8) (see Fig. [2]). Since the 
velocity profile should be stationary in the lab frame, the ions must be inhomogeneously 
spaced. Additional electrodes exerting a force F e (8) on the ions generate the necessary 
(de-)acceleration. The oscillating displacements of the ions around this equilibrium 
motion are phonons with velocities c (8) oc (v (8))~ L/ . Re gions with sufficiently large 
ion spacings and sufficiently low phonon velocities are supersonic; here phonons can only 



| After the submission of this paper the observation of the quantum Hawking effect has been claimed 
in [16] . But pair emission of and entanglement between Hawking particles has not been observed yet. 
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move in the direction of the ion flow and are trapped like light inside a black hole. We 
consider a system with a super- and a subsonic region. The border between these regions 
is analogous to a black hole horizon and will be shown to emit Hawking radiation. 

The plan of our paper is the following. First we will describe the system of ions on 
a ring and explain why it is expected to show the Hawking effect in Sec. [2] Then we 
briefly review previous works which are relevant for our analysis in Sec. [3j After this 
preparation of the reader we present our simulation results and discuss them in Sec. HI 
We are analyzing the Hawking effect in our system in two distinct ways: First we prove 
that the Hawking effect has a thermal spectrum by scattering pulses on the black hole 
horizon, second we are analyzing the emission of Hawking radiation from a black hole 
after its formation. In the latter situation we analyze the entanglement between pairs of 
Hawking particles, one inside and one outside the black hole. In this context we observe 
a transition from the quantum to the classical Hawking effect (spontaneous versus 
stimulated emission). Most importantly, this is suitable for experimental verification. 
In Sec. [5] we will discuss experimental setups which allow for the measurement of the 
analyzed physics. 



2. Ion Ring System 

In this section we are presenting our theoretical description of ions on a ring. In Sec. 12.11 
we start with the full Hamiltonian, explain our assumptions and approximations, and 
point to the equations underlying our simulations. We establish the connection with 
the general description of analog black holes in hydrodynamical system by considering 
the continuum limit of the system in Sec. 12.21 



2.1. Discrete Ion System 

In the following we are explaining the detailed setup of our proposal. The dynamics of 
N ions with mass m and charge e are described by the Hamiltonian 

i=l 1 i=l 

with the Coulomb potential V c and a local external potential V c (8). Instead of 
specifying V c we will impose an angular velocity profile v (9) by fixing the equilibrium 
positions. The required V e is then determined through the difference between the ion 
acceleration 9® (t) and the Coulomb force F? (9®, . . . , 9%). The Coulomb force and the 



external force are given in Appendix A 



To this aim we impose the classical equilibrium positions 

«?(') = + P) 

where g maps the normalized indices i/N G [0, 1] monotonically increasing onto the 
angles 9 e [0, 27r] and is periodically continued, g must be sufficiently smooth, i.e. three 
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times continuously different iable. The stationary angular velocity profile is 

*W = *^). (3) 

where T denotes the rotation time of the ions (see Fig. [2]) and g' is the derivative 
with respect to the argument of the function g. In a part of this paper we consider to 
dynamically create a black hole metric from a flat metric. To this aim, we decrease v m i n 
from the value f m ; n = 2it/T for homogeneously spaced ions in a Gaussian way with time 
constant r (see Appendix B). 



We choose a stationary velocity profile as depicted in [2j It is composed of a 
subsonic region with the constant angular ion velocity v (9) = v m i n in the angular 
range < 9 < crf min T and 2n — crt> min T < 9 < 2n. In the supersonic region 
o~v m i n T < 9 <2tx — av m i n T the constant angular ion velocity is 

n 1 _ 2o-t) min T 

v(e) = Vmax = - - 9 2 - . (4) 

i 1 — 2(7 

The two velocities f m j n and v max are naturally constrained by the rotation time of the 
ions v-m in < 2n/T < v max . The black hole horizon in our system is located close to 
#h = o-v min T. Due to the ring structure a second horizon, the white hole horizon exists 
at $h = 2tt — o~v m i n T. The transitions at the horizons between the subsonic and the 
supersonic regions contain 271 (black hole horizon) and 272 ions (white hole horizon). 
The exact expression for the velocity profile is given in the Appendix in Eq. (IB.ljl . 

We choose w m i n = 2ir x 0.83/T and that the small transition regions 271, 272 contain 
0.04Af and O.IA^ ions unless otherwise stated. If necessary, the black hole is dynamically 
created in the small time interval r = 0.05T (see Eq. ( IB. 211 ). 

In this paper we do not work with the full Hamiltonian (pQ). Instead, we treat small 
perturbations around the equilibrium motion 61 (t) = 9® (t) + 59i (t) and expand the 
Hamiltonian to second order in 89 i 

1 N 

i=l i^j 



with the time dependent force matrix J 7 = (/y) (see Appendix A ) and the canonical 
operators 8§i and L8pi/(2Tr) = —ihdsg t describing the phononic oscillations of the ions. 
This harmonic approximation is valid if the typical variation of the ion position is much 
smaller than the ion spacing 



^7 2tt / h A . 1\ 2tt A . 1\ 2tt 
" 2> " TV m^W ( <n> + 2) " (<"> + 2 J « N ' ^ 

Here we use the typical mode frequency Afa; rot where u TOt = 2n/T is the rotation 
frequency of the ions. If we tune the ion velocity to lie in the same order of magnitude as 
the sound velocity, ui rot gives the smallest mode frequency. We insert the experimental 
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parameters considered in Sec. with iV = 1000 here. In this paper we propose an 
experiment where the Hawking temperature Th and the initial temperature T are not 
much larger than hw ro t/kB, so we can assume (n) ~ 1. Then the above requirement is 
nicely fullfilled. 

The quasi-free quantum dynamics of this harmonic system (jSJ) are governed by the 
classical linear equations of motion for the first and second moments. Through Wick's 
theorem all higher order correlation functions are determined by these moments. The 
first moments are with the definition 




ie{0,...,N-i} 

i G {N, . . . , 2N — 1} 



(7) 



The second moments are grouped under the covariance matrix 

1 



r 



2h 



<{£<&}- 



(8) 



where {} + denotes an anticommutator. The equations governing the dynamics and 
determining the equilibrium states of the first and second moments are given in 



Appendix C 



A stability analysis of this system is described in Appendix D We find that 
exponential instabilities, though present in this system, are not important for the 
proposed experiment because it is performed during only one rotation period T. 



2.2. Continuum Limit 

In order to get some insight, we consider the limit of an infinite number of ions and 
formulate the analogy with the standard Hawking effect in this limit. 

We will first study the behavior of the dispersion relation of the ions at small 
wavenumbers (see Fig. [3] for the finite chain). We use the approximation given in [M] 
for small wavenumbers in the open Coulomb chain taking into account the long range 
Coulomb interactions. From 



2 / ka 
u{k) = ck x jl- - log I — 



(9) 



with the ion spacing a = L/N we get 

du{k) lo 1 
dk k 3 



and 



d 2 L0(k) 
dk 2 



l-flog(fe) 



1 + 



(10) 



3W1-IMt) 
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Figure 3. Phonemic group velocity c(k) in the flat subsonic region as a function of 
k for full Coulomb interactions (blue dashed line) and nearest-neighbor interactions 
only (green solid line). For nearest-neighbor interactions only the group velocity 
approaches a constant at small wavenumbers displaying the linear dispersion. For 
full Coulomb interactions a logarithmic divergence at small wavenumbers is observed 
for a finite system size. We use av m - m T — 2ir ■ 0.375, N = 1000, and e 2 /47reo = 
1.2591/(2AT) • mL 3 T~ 2 (see |AppendixB ) 



For k = In/L ■ n with 1 < n < n <C N this can be summarized as 

^„" and ^ I8 0. (12) 
dk k dk A 

For a given k, the phase and group velocity become identical in the continuum limit 
N/L — > oo as explicitly shown in Eq. (TTU]) . This is important for the definition of 
the Hawking temperature in the following, see Eq. ffTTj) . If group and phase velocity 
did not coincide, it would not be clear a priori how to determine the correct Hawking 
temperature. For example, in the scenario considered in [20J, the product of group and 
phase velocity enters the formula for the Hawking temperature. Since both velocities are 
identical in the continuum limit for our proposal, for the finite system we can determine 
the Hawking temperature at a given k from the group velocity. 

Now, we calculate the Lagrangian for the scalar field $ (0? (t) , t) = 56i (t) in the 
continuum limit. Here we can make an analogy with the standard Hawking effect as 
observed in [3]. Because of the equilibrium motion of the ions the kinetic energy K 
transforms according to 

2 I 2tt J ^\dt 

x ' i=l x 

= d 6 p -^> (flfcS + v (9) cV£) 2 , (13) 
Jo * 
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where we introduced the conformal factor p(9) = n(8) ■ mL 2 /(2ir) 2 with the density 
n (9) = N/(v (6) T). For an homogeneous system the potential energy V transforms to 



N N-l 
i,j=l k=0 

N-l N T n 

L z m 



k=0 n,m=l ^ ' 







^ d9^-(iD(-id e )$(9)f, (14) 

where D (8, k) = c{9) k + O (k 3 ) is the dispersion relation of the Coulomb chain. 
Assuming a slowly varying v (8), we can now formulate the Lagrangian for the ion 
system in the laboratory frame in the continuum limit 



£ = j d8^ (d t $ + v(8)de$>y - (iD (9, -id e ) ^ 



(15) 



This scalar field satisfying a linear dispersion relation at low wavenumbers with sound 
velocity 

for nearest-neighbor interactions, resulting from the actual form of the matrix given 
in Eq. ( 1A.4j) . is analogous to a massless scalar field in a black hole spacetime as first 
shown in [3]. Its quanta cannot escape a supersonic region with v (9) > c (9) like photons 
trapped inside a black hole. The horizon of this analog model is located at c (#h) — v ($h) 
with 0H ~ cri7 min T. Following [3], pairs of Hawking particles are emitted close to the 
black hole horizon with a black body distribution at the Hawking temperature 

k B T n _ k 2_ 2\ I = 3 g"(g- 1 (9)) , 

h 2tt 4nvd9 [V ° ' |H 4vrT g> {g~ l {9)) ^= e »' 1 ' 

The first equality defines the surface gravity k (see [21 E]), in terms of the Hawking 
temperature. The second equality is derived in reference [3], the third one results from 
the explicit forms for v (9) and c (9) for nearest-neighbor interactions. 

In the case of long-range Coulomb interactions, the calculation of the Hawking 
temperature T H (see Eq. fTl7|) ) is difficult because of the non-linear dispersion relation 
at low wavenumbers. In order to estimate the impact of this dispersion, let us introduce 
an effective wavenumber- dependent Hawking temperature. Because the ion spacings are 
inhomogeneous, we use a local density approximation. At each angle 9 we calculate the 
local density n(9) of the ions. We then extract a group velocity c(9) from the analogous 
homogeneous system with constant density and Ng = [n(9)} ions, where the square 
brackets denote rounding to the nearest integer. The group velocity is calculated from 
the dispersion relation at adjacent wavenumbers. With this method we calculate angle- 
and wavenumber-dependent group velocities c(9, k). By comparing c(9, k) with the ion 
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Figure 4. Dependence of the Hawking temperature Th on the discrete mode 
frequencies w(k) in the comoving frame at small frequencies. For nearest- neighbor 
interactions, the variation with the wavenumber is not visible, for long range Coulomb 
interactions variations are present. This is due to the nonlinear dispersion. Since 
these variations are not very large, we may approximately use a single Hawking 
temperature for the full system. We use av min T — 0.375, N — 1000. (red) Nearest- 
neighbor interactions with e 2 /47ren = 27vT^-; (black) Full Coulomb interactions with 
e 2 /47re = (see | Appendix B[ ). 



velocities v (8), we find for each wavenumber the black hole horizon #h- The derivative 
of c in Eq. (ITT)) is performed with respect to the ion number N$ instead of the angle 
9 so that only a single rounding procedure Ng = [n{6)] is needed. Thus, Eq. ( ITT)) is 
transformed into 



k B T R (k) 1 g"{g- l {9)) 



h 



2nT g'(g-H6)) 
i fig-HQ)) 

2nT g'(g-HO)) 



1 + 



0n 



l + N f 



9h 



g'(g-HO)) 2 

c (N dH + 1, k) 







dN, 



c{Ng,k) 



Na=Ng 



c(N eu ,k) 



g'(g-H0)Y 



The dependence of the Hawking temperature on the wavenumber is depicted in Fig. H] 
for the parameters used in the analysis in Sec. 14.11 For nearest-neighbor interactions, 
the variation with the wavenumber is not visible, for long range Coulomb interactions 
variations are present. This is due to the nonlinear dispersion relation. We make an 
arbitrary choice and use one of these temperatures at small wavenumbers to analyze 
the whole system at all wavenumbers. In Sec. 14.11 we will discuss how our results 
depend on this choice. Note that in the continuum limit the dispersion relation for full 
Coulomb interactions becomes linear which resolves the arbitrariness. Thus we find a 
single Hawking temperature for the full system based on the form of the velocity profile. 
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3. Review 



In this section we will give a short summary of previous works related to our proposal 
(see [25J for further details). In Sec. 13. II we are describing a basic derivation of Hawking 
radiation with quantum field theory in curved spacetime for systems with a strictly linear 
dispersion relation. Secondly, we present a theoretical method that tests whether the 
expected radiation has a thermal spectrum [2TJ [18] for systems with sublinear dispersion 
relations (see Sec. 13. 21) . Then we explain a proposal, originally aiming at Bose Einstein 
Condensates [6], on how to detect Hawking radiation in an experiment as quantum 
correlations which emerge between a supersonic and a subsonic region, i.e. between the 
inside and the outside of a black hole (see Sec. 13. 3j) . 



3.1. Hawking radiation and Mode Conversion 

In this section we will summarize the derivation of Hawking radiation for a massless 
scalar field with a linear dispersion relation D (k) = ck in a black hole analog spacetime, 
defined by a velocity profile v (9) (see Eq. ( IB. 31) and Eq. (IB. II) ) and a density profile 
p (9). This system is governed by the Lagrangian 



£ = J d9^p- (d t $ + v (9) d g &) 2 - c 2 d e & 



(19) 



We will in the following assume p (9) to be constant and work with the field $ = yfpQ. 
The correspondend classical field equation is 

[(d t + d e v (t, 9)) (d t + v {t, 9) d e ) - c 2 df\ * (t, 9) = 0. (20) 

The solutions of this equation determine the notion of excitation modes in the system 
and are of the form 

/■oo 

*M) = E / du(*i{t,0)+H.c) (21) 
± Jo 

with \^ (t, 9) = exp (±iut) ^ (9) for a stationary spacetime. u denotes the frequency 
in the lab frame, it is related by u; = (v ± c)k ■ L/(2ir) to the frequency ±ck in the 
comoving frame. The solutions ^ are normalized with respect to the Klein-Gordon 
inner product 

(^|^/> = — / d9 [** (d t + vd e ) ^ - * w i (d t + vd e ) %} . (22) 

For each frequency u> we can find four independent modes ^t/^t*- Complex 
conjugation relates modes with a positive /negative sign of the Klein-Gordon norm 
M w = r^w)j corresponding to positive/negative frequency modes or particles/anti- 
particles. Note that the frequency in the comoving frame ±ck determines the sign of 
M . The index ± denotes modes which are left-/rightmoving (up-/downstream) in the 
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comoving frame. Finally, the scalar field theory is quantized by expanding the field 
operator in the modes (we drop the summation index ± now) 

* = E + a M ( 23 ) 

Ul 

and postulating canonical bosonic commutation relations for the mode operators a w 

[a«,aw']=0, [a\ i ,a u i\=8 U j U >. (24) 

The vacuum of the system is defined to be annihilated by all a w 

a u \0) = 0, Vw. (25) 

We want to study the following time dependent situation: At initial times, called in, the 
black hole is not present, the vacuum state is denoted |in) with corresponding modes a™. 
Then a black hole is created. At final times, called out, the vacuum state of the system 
is | out) with corresponding modes a° ut . We are here interested in the time evolution of 
the system from initial to final times starting from the vacuum state of the system | in) . 
The result of this time evolution can be described by a Bogoliubov transformation on 
the classical modes or the quantum mode operators 

K ut = E (<w*5 + /u^ir) , (26) 

us 

c* = E - ( 27 ) 

with the Bogoliubov coefficients a^/ and Pu>u>'- For a static and continuous system 
the frequency u is conserved and the Bogoliubov coefficients are diagonal, e.g., = 
b~w,u>'Pu- The coefficients are related through ja^l 2 — \/3 u \ 2 = 1. The particle content of 
the state (out) in terms of the state |in) is related to the Bogoliubov coefficients through 

(in|iVr|m) = El/ 3 -'| 2 =l^| 2 - (28) 

Therefore, it is nonvanishing if the time evolution from initial to final times mixes 
positive and negative frequency modes. In particular, the production of Hawking 
particles can be understood as the evolution of initial positive frequency modes into final 
negative frequency modes. The Bogoliubov coefficients of the inverse transformation 
have the same modulus as the ones of the forward transformation 

*i n = E k^t - p<*s*r) . (29) 

^ = E (^C - ■ (3°) 



We will employ this relation to numerically determine the Bogoliubov coefficients for 
the discrete case of ions on a ring. 
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In the following we will briefly summarize the analytical calculation of the 
Bogoliubov coefficients for a continuous system with a strictly linear dispersion relation. 
One can decompose the field operator at late time into an upstream and a 
downstream ^ + part 

£ = + (31) 

The downstream part is not strongly affected by the black hole creation, hence the 
particle creation in this part is negligible /3+ = 0. This assumes that the mixing between 
upstream and downstream parts is small - otherwise we would have to include a grey- 
body factor, cf. So Hawking radiation is created by the upstream part only. It can 
further be decomposed into modes \&<J supBr inside the supersonic and ^/J sub inside the 
subsonic region. The decomposition of the field operators is finally 

*" = ^supcr + * S "ub, (32) 

POO 

Kupex = du) (Csuper^supere" 1 "* + H.C.) , (33) 

JO 

*«b = J™ duj (a-L^e-^ + H.c) . (34) 

We have to find the regular in field after the horizon formation ^~ in , decompose it in 
terms of the out modes 

Km (*) = a ^,sub (*) + /W,supc r (0) , (35) 

and read off the Bogoliubov coefficients a w and f3 w . So in summary, we calculated the 
time evolution of global upstream waves before the black hole creation into upstream 
waves after the black hole creation on both sides of the horizon. In the lab frame the 
resultant waves travel away from the horizon on both sides. 
The interesting Bogoliubov coefficient /3 W is found to be 

\^? = — Ur (36) 

-p {-&) - 1 

with the Hawking temperature Th (see Eq. f|T7|) ). We summarize these findings by 
stating that Hawking radiation emitted after the creation of a black hole has a thermal 
spectrum 

(m|iVr|m) = J2 7^71 : ■ (37) 

w exp [btfk 

3.2. Scattering of Pulses 

The derivation of Hawking radiation presented in the last Section l3~TTl suggests calculating 
the Bogoliubov coefficients from the time evolution of classical pulses scattering at the 
black hole horizon. It is numerically advantageous to calculate the scattering process 
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Figure 5. Schematical depiction of the dispersion relation for the scattering process 
of a final negative frequency pulse from the black hole horizon. Negative/positive 
frequency modes are depicted with solid/dashed lines. Upstream/downstream modes 
are depicted in red/green. The situation at late times is presented in (a), at 
intermediate times, when the pulse hits the horizon, in (b), and at early times in 
(c). Part (d) shows the situation outside of the black hole for small ion velocities, for 
which we observe the effect analogous to Bloch oscillations (see text). 



backwards in time. In the following we will focus on the situation of a sublinar dispersion 
relation in a continuum system as studied in [2T] . 

The calculation starts from a pulse with negative frequencies and small positive 
wavenumbers at late times. This pulse is moving upstream and is leftmoving in the 
comoving as well as in the lab frame. We then calculate its history in time: It 
approaches the horizon and is reflected by it. The reflected pulse mainly consists of 
two early-time upstream pulses, one positive and one negative frequency pulse, at high 
absolute wavenumbers. The positive/negative wavenumber pulse has negative/positive 
frequencies in the comoving frame, i.e., both are upstream pulses. But they have such 
a small group velocity - remember the sublinear dispersion relation (see [3]) - that they 
are dragged along with the ion background. So they are rightmoving in the lab frame, 
while being leftmoving in the comoving frame. These pulses can be determined due to 
frequency conservation for a static system in the lab frame. Note that the sign of the 
frequency in the comoving frame determines the sign of the Klein- Gordon norm (see 
Eq. (|22l) ) and thus the notion of particles versus anti-particles. Therefore, the relation 
between the late-time negative frequency and early-time positive frequency pulse is 
important for the particle producing Hawking effect. 
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Now, we will explain this mode conversion with the help of this frequency 
conservation and the dependence of the dispersion relation on the local ion velocity. 
The frequency in the lab frame is 



i.e., it is the Doppler shifted dispersion in the comoving frame Co> com = ±D(k). In 
Fig. [5] we depict the pulses at different times during the simulation (a-c) on the 
dispersion relation in the stationary lab frame together with the initial pulse frequency 
Uq. For comparison the interested reader can compare with Fig. E] depicting the 
pulses in real space. We are starting our discussion at late times (see Fig. G2a)) 
with a single pulse, upstream and with negative frequency. Its negative group velocity 
in the lab frame means that it is leftmoving in the lab frame. When the pulse is 
approaching the black hole horizon backwards in time the ion velocity at the pulse 
increases. In comparing Fig. E^a) and Fig. [5(b) one can observe the blueshifting of 
the pulse. When the pulse reaches the horizon, it is approximately at the minimum 
of the dispersion relation (see Fig. [5(b)). Now the mode conversion occurs, which 
is restricted by the frequency conservation in the stationary lab frame. Apart from 
the late-time pulse, three pulses are in agreement with frequency conservation. One 
pulse is moving downstream at low negative wavenumbers with negative frequencies. 
Thus, in the lab frame it is a fast rightmoving pulse. The two other solutions are 
at high absolute wavenumbers and are upstream, but slowly rightmoving in the lab 
frame. The positive /negative frequency contribution of the upstream pulse is located 
at positive /negative wavenumbers. Backwards in time all three pulses travel leftwards 
away from the horizon (see Fig [5(c)). Note that in [21] the downstream pulse is not 
observed in the actual dynamics. 

From a comparison between the Klein-Gordon norms of the late-time negative 
frequency pulse and the early-time positive frequency pulse the Bogoliubov coefficients 
1/3^1 2 can be extracted. In the literature a different approach is chosen, with the 
assumption that Hawking radiation is thermal, i.e. that the Bogoliubov coefficient 
is given by Eq. ( 1361) . the Klein-Gordon norm [22] of the early-time pulse is calculated 
from the late-time pulse. The Klein-Gordon norm of the early-time positive frequency 
pulse is 



^lab = vk— ± D(k), 





(39) 



It is compared with the prediction from the late-time pulse 




(40) 



Either the norms are compared in total or frequency-wise with the help of the expression 
dk = j^du, i.e. by division with the pulse velocities in the lab frame. One can compare 



Hawking Radiation on an Ion Ring in the Quantum Regime 



15 



the early-time positive frequency pulse 

M k (t < 0) 



Af k + = _ ' (41) 

^min Ck 



with the prediction from the late-time negative frequency pulse 



A? = ^ (i = 0) ' , (42) 

Vmin - C k pvn ( Ml 



<>X l> \ k B T 

and the prediction from the early-time negative frequency pulse 

~ N k (t < o) / hu k \ 

N k = exp -— - . (43) 

Vmin ~ C k V k BT J 

In [21] the thermal hypothesis is confirmed using both the integrated and the modewise 
comparison. 

In [18] a discretized hydrodynamic system is treated with the same method. On 
a lattice the dispersion relation is naturally sublinear. Since the lattice is moving, the 
frequency in the lab frame is not conserved anymore. It is shown numerically exact and 
with analytical approximations (WKB theory) that the mechanism of mode conversion 
described in [21] persists in this scenario. Even the Bogoliubov coefficients extracted 
from the comparison of early-time and late-time modes agree with the predictions for a 
continuum system with strictly linear dispersion relation (see Eq. fl36|) ). The existence 
of a downstream pulse is not reported in [18] , too. 

For a discrete and finite system with iV particles on a system of size L a finite 
number of wavenumbers exists 

k€ fJ N '-v*,J N, -^,...m (44) 



L 

with the renormalized ion number 

N' = N^l^. (45) 

v 

N' appears because the local ion spacing is not L/N, but a(8) = L/N' = L/N ■ 
v(9)T/(2tt) in the inhomogeneous system. This Brillouin zone is already visualized 
in Fig. |5j In this system a symmetry under combined translations in space and time 
survives, 

% ->• % + 1 (46) 
t->t + T/N (47) 

(see Eq. (j2j)). In a homogeneous part of the system this symmetry implies that the 
state of the system is invariant under the combined transformation 

k -> k + —N'n, (48) 

2tt . . 

u) -)> oj + —Nn (49) 
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for any integer n. 

We use this transformation to explain an effect described in [18] analogous to 
Bloch oscillations. The two high wavenumber solutions of the frequency conservation 
condition can lie outside the Brillouin zone. Then we find the solutions by looking at the 
frequencies uo ± 2nN/T, depicted in Fig. [5(d). Thus, the early-time pulses are instead 
located on the downstream branch of the frequency condition. The positive/negative 
frequency pulse is located at high positive/negative wavenumbers (see Fig. [5](d)). 

3. 3. Correlations 

The Hawking effect produces pairs of particles propagating away from the horizon in 
opposite directions, one outside the black hole (the Hawking particle), and one inside the 
black hole (its infalling partner). Taking advantage of this, it has been proposed to detect 
Hawking radiation via correlation measurements (see [6J), which reveal the entanglement 
between the two Hawking partners. In [6J the density-density cross-correlation of 
phonons on both sides of the horizon produced through the Hawking process in a sonic 
black hole built up with a Bose-Einstein condensate (BEC) is studied [5]. In a quasi one- 
dimensional weakly interacting BEC in the hydrodynamical approximation the sound 
propagation is described by 

C = J d9^ + v (9) d e $>) 2 - c 2 d e § 2 ] . (50) 

This is manifestly the same as Eq. ( 1T91) apart from the actual form of the conformal 
factor /C, which if constant does not affect the dynamics appart from a rescaling of 
the fields. So the propagation of phonons in a BEC, i.e. phase/density excitations, 
and the propagation of phonons in a ring of ions, i.e. displacements of ions from their 
equilibrium position, are described by the same physics, and therefore share similar 
behaviors. Thus, the analysis performed in [6] can be translated into the context of ion 
rings. 

The general equal-time two-point correlator (\& (9)^f (9')) in 1+1 dimensional space- 
times is given by [27] 

<* (0)* (0')> = - lim A log [(17(0, t) - U(9',t')) (V(9, t) - V(9', If))] . (51) 

This is the stationary correlator of a massless minimally coupled quantum scalar field in 
1 spatial dimension. It exhibits a characteristic structure with the logarithmic distance 
between the two points as the difference between the null coordinates U and V [27]. The 
meaning of the null coordinates is the following: the mode solution can always be written 
in the form ~ e ~ lujU f or upstream modes and ~ e~ tulV for downstream modes. The modes 
propagating upstream (downstream) move along space-time trajectories that keep the 
V (U) coordinates constant. For example, if the system is at rest and homogeneous, 
the analog spacetime is the Minkowskii (flat) one, and the Kruskal coordinates reduces 
to the familiar U = t + 9/c, V = t — 9/c [29], that define the standard light- (sound-) 
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cones. If the system is more complicated (i.e. non-homogeneous or moving), the mode 
propagation is different, and the associated null coordinates will display distortion. 

Eq. ( 15T1) applies to all conformally invariant theories in 1+1 dimensional spacetimes 
(which are always conformally flat). In the present case, the theory we are dealing with 
is not conformally invariant for the presence of the conformal factor K, in the action. 
Nevertheless, assuming that K, varies smoothly over the system, the correlator can be 
approximated using Eq. fl5Tj) also in the present case by 



(He)He')) = - iim A log [(u(9, t) - u{e\ t')) (v(e, t) - v(9', t'))) , 

(52) 

where the 1/ yK,(0)1C{0') term follows from a rescaling of the field, and terms containing 
derivatives of /C are neglected. 

Let us move to the evaluation of this correlator in the presence of Hawking radiation. 
Since the modes responsible for the emission of Hawking radiation are the upstream 
modes we will focus on the upstream sector of Eq. ( 15T1) only. The equal-time 
two-point correlator in the presence of a black hole horizon is given by [27J 



<*-(*)*-(«-)> = - lim A_l_ log lm t) _ v(f/< i0) . (53) 

The downstream part of the correlator remains unaffected even in the presence of a 
horizon, i.e. V still reads t — 6/c. The only modes which get distorted by the presence 
of the horizon are the upstream ones which become far away from the horizon 



U(9, t) m/out = ± e -^+«5j ) (54) 

for the interior (+) and the exterior (-) region of the black hole. They suffer the 
typical exponential distortion of the upstream modes due to the presence of a black 
hole horizon [27] 13] (discarding transients). The exponential distortion follows from the 
wave equation for a linearized velocity profile at the horizon, k is the surface gravity 
on the sonic horizon which is proportional to the Hawking temperature (see Eq. (fl7"|) ). 
It is worth emphasizing that the form of the modes (1541) is universal for any black hole 
horizon formation, independent of the details of its formation. This is the origin for the 
universal behavior of Hawking radiation. 

The momentum-momentum correlator (Sp~(6)Sp~(8')) can be obtained from Eq. 
(|53|) . The conformal factor /C must be replaced with the conformal factor from the 
Lagrangian of ions on a ring (see Eq. (1191) ). With the relation 

5p~ = Am(<9 t + vd e )&~ (55) 

27T 
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one gets 



{8p-(6)5p-(6>)) 



hm 



16tt [c (0) - v {9)} [c (0') - v {6')} 



c{6)c{9>) 



K 



n(9)n(6') cosh 2 



c(e)-v(p) c(fi')-v{p') 



■ (56) 



This correlator has the typical form associated to the Hawking effect. The correlations 
are scaling with the square of the Hawking temperature. For and 9' on opposite sides 
of the sonic horizon, c (9) — v (9) and c (9') — v (9') have opposite sign. Therefore, the 
momentum-momentum cross-correlations are negative. They exhibit a peculiar peak 
along a straight line for 9 = ^y"}^,^ ' oc 9'. These cross-correlations correspond to the 
two entangled Hawking particles propagating in opposite directions as they move apart 
from the horizon. 



4. Simulations for Ion Rings 

We are now returning to the discussion of the discrete ion chain. In this section we are 
presenting the results of our simulations and are comparing them with the predictions 
and expectations from Sec. |3j We are pursuing two routes of simulations: First we 
simulate the scattering of pulses on the black hole horizon in Sec. 14.11 From the result we 
can deduce the Bogoliubov coefficients and theoretically confirm the thermal hypothesis 
(see Sec. 13. 2j) . A second series of simulations presented in Sec. 14.21 is analyzing the 
emergence of correlations between the inside and the outside of a black hole after its 
creation (see Sec. I3.3p . These correlations demonstrate the pair creation mechanism of 
Hawking radiation and are closely related to the emergence of entanglement between 
the inside and the outside of a black hole. In contrast to the scattering analysis, these 
simulations also act as direct proposal for an experiment as further discussed in Sec. EJ 



4-1- Scattering of Pulses 

Our numerical results for the propagation of a final pulse backwards in time in the 
discrete system of phonons on an ion ring are presented in the following. We first 
introduce the quantities necessary for this analysis. If the phononic excitations are 
localized in the flat subsonic region with constant ion velocity t> m j n , the excitations 89i(t) 
and 59i(t) can be expressed as modes 56k(t) and b~9k{t) with wavenumber k, where the 
dot represents the time derivative in the comoving frame. Due to the finite system size 
only discrete wavenumbers appear (see Eqs. (jUl) and ( )45i) ) and thus the dispersion 
relation is also discrete. The positive and negative frequency part of these excitations 
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Figure 6. Ion displacement \S6i(t) / (2tt)\ during propagation backwards in time 
starting from the final wavefunction in Eq. (1611) with s = 5. One late-time pulse 
and three (two of which overlapp) early-time pulses are present (see text). Two early- 
time pulses have large absolute momenta (blue dashed lines) , one early-time pulse has 
small momentum (green straight lines). We use av m i n T = 0.375 and TV = 1000. For (a) 
and (b) e 2 /4^e = i ff il #, for (c) and (d) e 2 /4^o = (see | Appendix B) . 

(a) and (c) use lab frame angles 9, (b) and (d) use ion indices i (see text for description 
of scattering process. 

are defined by 

(*) = \ ( 50 k(t) ± i59 k (t)/u k } , (57) 
<^fc (*) = \ (Mk(t) T ioo k 56 k {t)) . (58) 

These relations follow from the spatial behavior 58 (6) ~ exp (ik^-L) in regions of 
constant ion velocity v (9) = v [21]. The analysis of the particle production requires us 
to use the Klein-Gordon norm for these modes. In our special case the Klein-Gordon 
norm defined in Eq. (|22j) becomes up to a constant 



fe,± 



(59) 
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Figure 7. Klein-Gordon norm A4 of the phonon pulses for propagation backwards 
in time starting from the final wavefunction in Eq. (|6f I) with s = 5. With a 
comparison of these pulses we test the thermal hypothesis. We use (Tw m i n T = 2tt- 0.375, 
and N = 1000 (see |Appendix B| . The final negative frequency pulses are depicted 
in green, the initial negative frequency pulses in blue (dash-dotted line), and the 
initial positive frequency pulses in red (dashed line). We depict the discrete norm 
distributions/dispersion relation as a continuous curve here, (a) Final (t — 0) and 
initial (t = -0.67T) norm distribution of S9 k (t); e 2 /4Tre = 1.2591/(2JV) ■ mL 3 T~ 2 . 
(b) Final (t = 0) and initial (f = — 0.51T) norm distribution of 56f.(t) showing Bloch 
oscillation; e 2 /47re = 2.0004/(2iV) • mL 3 T~ 2 pQ. 



with 

■A/? = Mi* 50 t - Sd k*^k- (60) 

We are now describing the numerical calculation of the Bogoliubov coefficients with the 
method presented in Sec. 13.21 In summary, we are calculating the history of a negative 
frequency pulse on the upstream branch of the dispersion relation that travels away 
from the horizon. Back in time it scatters off the horizon and originates from several 
pulses. The early-time positive frequency pulse and the late-time negative frequency 
pulse are related through the Bogoliubov coefficients. 
We start from the final pulses 

56° k (0) = fc-e"(^) 2 , s = l,...,20, (61) 
centered at different wavenumbers k to test different frequency ranges. We calculate 



its history with Newton's equations of motion given in Appendix C (see Eq. (jC.ip ) by 
using an iterative differential equation solver. From frequency conservation in the lab 
frame we expect three pulses on the upstream branch of the dispersion relation [2lT [18] 
(see Sec. 13. 2p . 

Before coming to the simulation results, we will discuss the parameter regime used 
in this section. The ratio between the typical pulse frequency u in the comoving 
frame and the Hawking temperature is of order fiw/ksTx ~ 10 (see also Fig. H]). 
In [30] the deviations from thermality have been examined based on the quantity 
w max = max fc>0 (D(k) — vkLjln). Here o; max takes the following values: uj max = 290/T 
for nearest-neighbor interactions only (Fig. IHk), Umax = 610/T for full Coulomb 
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interactions (Fig. [Hb). The surface gravity for the former case is k = 65/T, for the 
latter k = 82/T (see Fig. @J, thus the ratios are uj max /n = 4.5 and u) max /K = 7.4. 
This regime is identified as the regime of small deviations from thermality in [30]. This 
prediction applies to systems with nearest-neighbor interactions only. One of the main 
results of our paper is to consider also the long range Coulomb interactions for which 
the deviations from a linear dispersion relation at small wavenumbers are significant in 
a finite system. 

Our results agree mainly with those of references [2T| [18] as shown in Fig. [7J 
First, we observe quantitative devisions, whose order of magnitude agrees with the 
uncertainties in calculating the Hawking temperature. Second, we find a previously 
undescribed downstream pulse at small negative frequencies. But before discussing the 
detailed analysis of the scattering process, we will repeat the description of the mode 
conversion from Sec. 13.21 but this time illustrated in real space with the simulation 
results from the special case of the ion ring (see Fig. [6]). 

We begin with the normal scenario of sufficiently large ion velocities (see Figs. 
[6(a)- (b) and [7(a)). We observe all three early-time pulses which are in agreement with 
frequency conservation. Two pulses are located on the upstream branch at high absolute 
wavenumbers (red and blue lines in Fig. [7(a)). But at these wavenumbers they have 
such a small group velocity in the comoving frame that they are moving rightwards in 
the lab frame (blue dashed lines in Figs. [6(a)- (b)). A third pulse is located on the 
downstream branch of the dispersion relation at small negative wavenumbers (blue line 
in Fig. [7(a)). Its group velocity is large both in the comoving frame and in the lab 
frame. This time evolution is depicted in Fig. [6(a) in the lab frame and in Fig. [6(b) 
in the comoving frame. The two upstream and the one downstream pulse can clearly 
be identified in the comoving frame. In the lab frame the downstream pulse is quickly 
moving rightwards, the two upstream pulses are slowly moving rightwards. 

A different situation arises for sufficiently small ion velocities (see Figs. [6(c)- (d) 
and [7(b)). In this case, all three early-time pulses are moving downstream. This time 
evolution is depicted in Fig. [6(c) in the lab frame and in Fig. [6(d) in the comoving 
frame. In the comoving frame we can now observe the three downstream pulses, one 
fast pulse which is the one at low wavenumbers (green lines) and the two interesting 
slow pulses (red dashed lines) which are still located at high absolute wavenumbers. In 
the lab frame the two pulses at high absolute wavenumbers are travelling faster than in 
the normal case. 

The spectral analysis in Fig. [7J confirms these explanations. We find the two 
previously mentioned high wavenumber pulses. For small ion velocities the effect 
analogous to Bloch oscillations is observed. Then the main incoming pulses have 
wavenumbers with opposite signs. In addition we observe the weak downstream pulse 
with small negative wavenumbers that is not described in the literature. We find that 
its magnitude depends on the strongly on how the fringes of the pulse in outside the flat 
subsonic region are treated. Therefore we suggest that the presence of this pulse could 
be a consequence of the finite initial and final excitation probability of ions outside the 
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Figure 8. Comparison of spectrally resolved Klein-Gordon norms AF2 (late-time 
negative frequency pulse) at t = with A/jf (early-time positive frequency pulse) 
and AfjT (early-time negative frequency pulse) (see Sec. 13. 2|) as a function of the lab 
frame frequencies at t — —2T after propagation backwards in time starting from the 
final wavefunction in Eq. (|6Tj) with s — 5. We use av m i n T = 0.375, N = 1000. (a) 
Nearest- neighbor interactions at t = —2T with e 2 /47reo = 3 ) Full Coulomb 

interactions at t = — 0.75T with e 2 /47reo = 2^ 2 ^4- (see | Appendix B[ ). 

flat subsonic region. 

We further compare the Klein-Gordon norm of the positive frequency early-time 
pulse with the prediction for thermal radiation (see Sec. 13.21 and [2TJEEE])- We calculate 
the Hawking temperature according to Eq. (ITS]) . If only nearest-neighbor Coulomb 
interactions are considered, the relative difference between these norms is lower than 
e = 0.01 for N = 1000 ions. For the long range Coulomb interactions a conservative 
estimate yields e < 0.2. The latter result holds for any arbitrary choice for the Hawking 
temperature at the wavenumbers k = 2ir/L4 . . .6. The bound on e agrees with the 
differences between analysis based on the group and the phase velocity of the phonons 
(see Sec. Oand Eq. flT^). 

We also perform the analogous spectrally dissolved comparison (see Sec. 13.21 and 
[21]), comparing the early-time positive frequency pulse with the prediction based on 
the thermal hypothesis. Only discrete frequencies appear in the system because only 
discrete wavenumbers are present. Therefore, it is in general not possible to exactly 
match the frequencies of the pulses. Especially, we cannot perform a thermal fit. Thus, 
we look at the two frequencies in the early-time pulses closest to the frequencies in 
the late-time pulse. We compare the early-time positive frequency pulse J\f£ with the 
predictions for it based on the late-time negative frequency pulse J\fj? and based on the 
early-time negative frequency pulse AffT (see PT]) -( 143|) ). The result of this analysis is 
shown in Fig. Ufa) for nearest-neighbor interactions and in Fig. Mjo) for full Coulomb 
interactions. For this analysis we choose the Hawking temperature at k = 2n/L ■ 5 
(see Eq. (JTHjl ). but similar results are obtained for the adjacent wavenumbers). This 
analysis again confirms the thermal hypothesis to the extend possible. The accuracy of 
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Figure 9. Momentum- momentum correlations Cij(t) at time t = 0.5T with real space 
lab frame positions. Starting from homogeneously spaced ions with temperature T = 
a black hole region is created in the small time interval r = 0.05T. The dashed lines 
point to the signature corresponding to the emission of pairs of Hawking phonons. We 
consider N = 1000 ions, and av mm T = 2ir ■ 0.25. (a) Nearest-neighbor interactions 
with e 2 /47re = ^ (b) Full Coulomb interactions with e 2 /47re = ° 2 2 ^ 53 m T L " [J] 
(see |Appendix B| 



our analysis is restricted by the discreteness of the system and the nonlinearity of the 
dispersion relation at small wavenumbers. 

We summarize the findings of this section: Hawking radiation with a thermal 
spectrum is emitted from a black hole horizon on an ion ring, even for a finite system 
with a logarithmically diverging group velocity at low wavenumbers due to long range 
interactions (see Fig. [3]). 

4-2. Correlations 

For an experimental proof of Hawking radiation on ion rings, we propose to observe the 
emission of Hawking radiation following the creation of a black hole horizon (see Sec. [5]). 
We propose to measure the emitted phonons or the emerging correlations between the 
subsonic and the supersonic region (see [6]). The latter is discussed in Sec. 13.31 in this 
section we are presenting simulation results for ion rings on the emerging correlations. 
We compare our results to the analytical findings for a continuum system with linear 
dispersion relation derived in Eq. ( |56|) [Sj. In the quantum regime the emergent cross- 
correlations display the generation of entanglement. We study its properties, especially 
to analyze the crossover from the quantum to the classical Hawking effect. 

4-2.1. Discrete System We propose the experiment to start from the ground/ thermal 
state of the excitation around homogeneously spaced ions at rest with temperature 
To. Then the system is accelerated with a constant force which does not change the 
quantum state of the system defined relative to the equilibrium positions. Subsequently, 
a supersonic region is created in the small time interval r. This is done by reducing 
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the subsonic fluid velocity v m [ n in a Gaussian way, while leaving the average rotation 
velocity constant (see |Appendix B|and Eq. (]B.2j) ). We reduce excitations created at the 



white hole horizon at 6/2ir »s 1 — av m i n T by a wider transition region at this horizon. In 
an experiment the magnitude and velocity of these excitations can be watched through 
careful measurements. 

We are analyzing the momentum-momentum correlations 

C tj = (dp.dpj} -T/(hm). (62) 

In the continuum limit these momenta correspond to the time derivative of the scalar 
field $ 

5p i ^(d t + v(9)dg)$(9,t). (63) 

The dynamics of the momentum-momentum correlations are given by the dynamics 
of the covariance matrix (see Eq. (1C.3j) ). In Sec. 15.11 we are explaining how these 
correlations can be measured in an experiment. 

Fig. [9] shows the simulation results for CV,. It displays a fixed time after the black 
hole formation starting from the ground state, i.e., T = 0. Correlations between the 
inside and the outside of the black hole are created close to the black hole horizon and are 
moving away from it as expected. As stated earlier these correlations correspond to the 
pairs of Hawking particles. We interpret their pure existence as a signature for Hawking 
radiation. Fig. E(a) shows the simulation for interactions between neighboring ions 
only, Fig. [9(b) for long range Coulomb interactions. The correlations behave similarly 
in both cases (see discussion below). They can still be observed for initial temperatures 
two orders of magnitude above the Hawking temperature [6] as shown in Fig. [10J In 
Sec. 14.2.31 we find that the Hawking effect is still quantum at such initial temperatures 
(see Fig. [T4"]) . In contrast to the entanglement, the cross- correlations actually remain 
present at arbitrarily large initial temperatures [26] (note T oc T for k^T^ 3> fojj^ see 
Eq. (JOl). 

The most significant correlation signal (i), i.e. the line of negative cross-correlations, 
demonstrates a basic property of Hawking radiation, it corresponds to two upstream 
phonons, one inside and one outside the black hole. An additional correlation feature 
(ii) is fully inside the black hole, corresponding to a pair of upstream and downstream 
phonons inside the black hole. The features (i) and (ii) have already been reported for 
analog black holes in a BEC (see J6]). The propagation velocity of these correlations 
c{6) ± v{6) depends on the group velocity, i.e. the dispersion relation, and the ion 
velocity. The angle of the cross-correlation signature (i) is determined by the ratio of 
the phonon velocities inside and outside of the black hole (c(#) — v(9)) / (c(9') — v{9')). 
The dashed lines in Fig. [9] show the predictions for the direction of the cross-correlation 
signal (see Sec. 13 .31) . 

In contrast to Fig. |9(a), which is very similar to black hole analogues in a BEC 
[6], the correlation plot for full Coulomb interactions in Fig. |9(b) displays a more 
complicated structure. The cross-correlation signal, is broader and there are additional 
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Figure 10. Momentum-momentum correlations Cij(t) for initial temperature Tq = 
102Th at time t = 0.5T with real space lab frame positions. We consider N = 1000 ions, 




lines of oscillating correlations. With higher resolution and/or for larger Hawking 
temperatures these effects can also be observed for nearest-neighbor interactions. We 
attribute these changes to the more complicated and sublinear dispersion relation. 

Even though the precise determination of the correlations for a nonlinear dispersion 
relation is a rather complicated issue (see, e.g., [31]), we may obtain a qualitative 
understanding by means of the following simple picture. Since group and phase velocity 
nearly coincide at small wavenumbers, the correlations are created near the horizon as 
a nice 1/ cosh 2 -pulse as in the case of a linear dispersion relation (see Eq. (1561) ). 

However, as this pulse propagates away from the horizon, the non-linear dispersion 
relation deforms it. For a fixed t', 9', the two-point function (5p(t, 9)5p(tf, 9')) obeys the 
same wave equation as 5p(t, 9) itself. In a homogeneous region of a stationary spacetime 
this is the same as for the field \l/ (see Eq. (120]) ). 



Since modes with larger k propagate slower than those with smaller k, a pulse with 
an initial 1 / cosh 2 -shape will be deformed during the time evolution similar to Fig. [TTJ 
The main pulse (global maximum) becomes broader and oscillations develop, trailing 
the main pulse (local minima and maxima), which are caused by the slower modes with 
short wavelengths. This deformation applies to the outgoing Hawking radiation and the 
infalling partners in the same way as both are moving upstream. In a finite system the 
deviations from a linear dispersion are significantly larger for full Coulomb interactions 
than for nearest-neighbor interactions. As a result, this simple picture explains the 
difference between Figs. [9](a) and[9](b). 

4-2.2. Comparison with Continuum System In Sec. 13.31 we have calculated and 
analyzed the cross-correlation signal for a continuum system with strictly linear 




(64) 
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Figure 11. Schematic depiction of the pulse propagation on an ion ring. The left 
graph depicts an initial Gaussian pulse profile, the right graph the dispersed pulse 
after time t. Trails of oscillations are following the dispersed pulse. 



dispersion relation. We observe these correlations also for our discrete system with 
sublinear dispersion relation as shown in Fig. |9j In this Section we are quantitatively 
comparing the simulation results for ion rings with the analytical results for a continuous 
system, focussing on the peak magnitude of the cross-correlations (see Eq. (|56|) ). 

We perform this comparison for varying Hawking temperatures. It is tuned by 
changing the widths of the horizon region 71 (see Eq. ( IB. II) ) keeping constant the other 
parameters. The comparison is shown for nearest-neighbor interactions in Fig. fT2T a) 
and for full Coulomb interactions in Fig. fT2T b). where the dots represent the simulated 
peak magnitude for the ion ring and the curve represents the analytic peak magnitude 
for the continuum system with a linear dispersion relation. 

For nearest-neighbor interactions, the results agree very well at small Hawking 
temperatures; for full Coulomb interactions, the agreement is quite good at moderate 
Hawking temperatures. The deviation at large Hawking temperatures is caused by the 
discreteness of the system. For large Hawking temperatures the horizon region, which 
almost completely determines the properties of Hawking radiation, is small and contains 
only a few particles. Then the discreteness of the system becomes relevant. In the limit 
of very large Hawking temperatures the peak height should mainly depend on the lattice 
spacing [26J. 

4-2.3. Entanglement Generation In this section we are discussing the creation of 
entanglement between the inside and the outside of a black hole following its creation. 
Entanglement is unique to quantum processes. Therefore, the existence of entanglement 
proves that one can observe the quantum version of the Hawking effect with our 
proposal. In contrast, correlations between the inside and the outside of the black 
hole are present both for initial thermal states T S> 2h in the classical regime and for 
initial quantum states T = [26J. So we analyze the crossover between classical and 
quantum Hawking radiation (stimulated versus spontaneous emission). Furthermore, 
our analysis emphasizes the importance we assigned to the cross-correlations. The 
emerging entanglement can be measured on two routes, either by measuring the 
covariance matrix through a measurement of correlation in the ion displacements (see 
Sec. 15. ip or by swapping the entanglement from the motional to the internal degrees of 
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2jc k B T H /h-T ~ ~ 2jtk B T H /hT 



Figure 12. Peak Height Cij(t) of the cross-correlation signature in the momentum- 
momentum correlations (see Fig. [5]) as a function of the Hawking temperature Tjj. The 
dots are simulation results for a discrete system with nearest-neighbor interactions only, 
N = 1000 ions, av m - m T = 27T-0.25, and r = 0.05T. 71 is varied to get different Hawking 
temperatures Th- h = IttK is Planck's constant, (a) Nearest- neighbor interactions 
with e 2 /47re = ^r 1 ^ and 0.003 < 71 < 0.1; (b) Full Coulomb interactions with 
e 2 /47re = 2 2 ^ 53 m ^ , and 0.005 < 71 < 0.1 (see |AppcndixB[ 



freedom of the ions 

We are now briefly introducing the relevant entanglement measures before we 
present numerical results for the ion system. The covariance matrix T (see Eq. (jSj)) 



can be calculated (see Appendix C ) and measured (see Sec. [5]) for the ion system in 
harmonic approximation (see Eq. (JSJ). It gives access to two entanglement measures: 
the entropy of entanglement and the logarithmic negativity. 

We consider a system with density matrix p divided into subsystems A and B. The 
entropy of entanglement is defined for a bipartite pure state. It is the Von-Neumann 



entropy of the reduced density matrix p^ of one subsystem A 

S(p) = -tr (p A log 2 p A ) • (65) 

The state p is a product state for S(p) = 0. As an entanglement measure S 
cannot increase under local operations and classical communications (LOCC). For 
pure Gaussian states S(p) can efficiently be calculated from the covariance matrix Y 
[3H [351 ES]. It is given by 

s (n = E ( A « A « - ^ - x ) ( A « - x )) ( 66 ) 

n=l 

with the symplectic eigenvalues X n , n = 1, . . . , N of the covariance matrix T. They are 
the eigenvalues of iaT with the symplectic matrix 



a 




(67) 
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Figure 13. (a) Time evolution of the entropy of entanglement S(T) for different 
Hawking temperatures Tjj with full Coulomb interactions (same parameters as Fig. 
|T2l^b) ) . (b) Rate of increase of the entropy of entanglement dS(T)/dt as a function 
of the Hawking temperature Xh- Full Coulomb interactions (red dashed line, same 
parameters as Fig. [T"2l b)). and nearest-neighbor interactions (black straight line, same 
parameters as Fig. IT2l a')'). The linear increase in entanglement is caused by the 
constant emission of entangled Hawking phonons. 



which exchanges position and momentum of each mode. 

We can also use the logarithmic negativity, which is an entanglement monotone, 
i.e., it does not decrease under LOCC [31]. In contrast to the entropy of entanglement, 
the logarithmic negativity can be calculated efficiently from the covariance matrix even 
for mixed states. It is defined as the logarithm of the trace-norm of the partial transpose 
of the density matrix 

N(p) = log 2 ||p T i (68) 

with ||M|| = yMXM. The covariance matrix r Ta of the partial transpose of p, i.e., p Ta , 
follows from the covariance matrix r of p by multiplying with —1 all matrix entries of 
T which contain exactly one momentum operator of subsystem A. Let k = 1, . . . , N 
denote the symplectic spectrum of T Ta . Then the logarithmic negativity is 

TV 

N(r) = - £ log 2 (2A n ) . (69) 

n=l 
2A„<1 

A vanishing logarithmic negativity does not mean that the system is not entangled, but 
it means that such systems cannot be purified to maximally entangled states. 

Now, we analyze the entanglement properties for black holes on ion rings. Using 
the entropy of entanglement we study how the entanglement is generated between the 
inside and the outside of a black hole starting from the ground state (T = 0) at initial 
times. The entropy of a black hole is more fundamental than the logarithmic negativity, 
but it is only defined for pure states. Thus, we have to compare the whole supersonic 
with the whole subsonic region. On a ring this means that one cannot determine at 
which horizon the entanglement is created. 
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t/T 



Figure 14. Time evolution of the logarithmic negativity N of a region with 0.2N ions 
(N = 1000) adjacent to the horizon for different initial temperatures To. The linear 
increase in entanglement is caused by the constant emission of entangled Hawking 
phonons. Full Coulomb interactions considered (same parameters as Fig. EJb)). 

In Fig. [l~3f a) we find a linear increase of the entropy of entanglement in time after 
an initial period (t ~ 0.2T in our case). This linear increase corresponds to the constant 
emission of Hawking radiation from a black hole. We have plotted the rate of this 
increase in Fig. fT3T b). We find for nearest-neighbor interactions a linear dependence 
of the entropy of entanglement on the Hawking temperature of the black hole. For full 
Coulomb interactions we find a similar behavior. 

The logarithmic negativity gives access to the entanglement developing for systems 
with finite initial temperatures T > 0. The entanglement between two small regions 
adjacent to the black hole is presented in Fig. HU The initial entanglement between 
the regions depends on the system temperature: For quantum systems at T = 
entanglement is present, for classical systems To — > oo it cannot be detected. The 
entanglement is increasing linearly in time for sufficiently small To. This behavior is 
already described above for the entropy of entanglement (see Fig. [TST a)). 

However, the logarithmic negativity allows further observations. For initial 
temperatures T more than two orders of magnitude above the Hawking temperature 
Th, no entanglement generation is visible in the logarithmic negativity. In contrast, 
the cross-correlation signal remains present for arbitrarily high initial temperatures. 
We observe here the transition from the quantum to the classical Hawking effect. We 
can conclude that for the initial temperatures T < Tq m IOOTh (N = 1000) the 
observed Hawking effect is quantum, whereas one would naively expect this transition 
at Tq w Th- Tq/Th is increasing with the number of ions, taking, for example, N = 100 
ions and e 2 /4vre = 0.6/(2iV) • (mL 3 )/T 2 ) we find T C « 24T H . T C /T H is approximately 
proportional to the largest mode frequency of the system, which is sublinear in in N. 

This might be understood with the following argument: As described in Sec. 13.21 
Hawking radiation emerges from large wavenumbers before being emitted at small 
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wavenumbers. The frequency related to these large wavenumbers is about iV-times 
(number of ions) higher than the smallest frequency in the system. Along this line, the 
Hawking effect remains quantum for initial temperatures comparable to these highest 
frequencies of the system at large wavenumbers. 

In Fig. [HI a saturation in the logarithmic negativity at later times is observed. 
This is in agreement with the fact that Hawking radiation is constantly emitted. After 
fully penetrating the small regions adjacent to the horizon, no additional entanglement 
can develop between these regions. 

5. Experimental Realization 

So far we have studied the appearance of Hawking radiation on ion rings. We have 
demonstrated that the emitted radiation has a thermal spectrum with a geometrically 
justified Hawking temperature (see Sec. 14. ip . and that it is created in pairs (see Sec. 
14. 2p . We complete our analysis in this section by describing an experimental setup which 
will allow for the observation of Hawking radiation on an ion ring. 

First, we present a suitable parameter regime for an experiment. The main 
condition on the experimental parameters is that the ion velocity must be approximately 
equal to the phonon velocity. This condition leads to the requirement e 2 /4ireo ~ 
0.25/ (2N)mL 3 /T 2 . For N = 1000 singly charged 9 Be ions with an average spacing 
of L/N = 2\im the rotation frequency of the ions would be u rot = 2n x 120kHz. oj TOt 
represents the smallest mode frequency of the system. For v min = (2n x 0.83)/T, 71 = 
0.02, and av m \ n T = 2n ■ 0.25 the Hawking temperature is h-aT^/H ~ 5/T pa 2tt x 95kHz. 
If the initial temperature is two orders of magnitude above the Hawking temperature 
To < IOOTh, we show explicitly in Sec. 14.2.11 that the cross-correlation signature of 
Hawking radiation remains present (see Fig. [10]) and we find in Sec. 14.2.31 that the 
Hawking radiation remains a quantum effect. Thus, it is not necessary to perform 
ground state cooling of all vibrational modes of the ions in an experiment. 

Note that it has been demonstrated long ago how to trap ions in quadrupole ring 
traps [22] and measure their arrangement [23]. The ideas of these experiments can be 
combined with modern cooling techniques applied to ions in linear Paul traps or in 
surface traps [371 EH] ■ Thus, the proposed experiment will allow to measure signatures 
of Hawking radiation for acoustic black holes with parameters and temperatures which 
can be reached in current experiments. 

The general idea of the actual measurement process is the following. The Hawking 
effect is encoded in the motional degrees of freedom of the ions, which are described 
in this paper with the ion displacements S9i and can be viewed as phonic modes. A 
different degree of freedom for ions is their internal state, here we address two hyperfine 
states of the ions. Lasers couple the motional degrees of freedom to the two relevant 
internal states. In this way, the information on the Hawking effect is transfered to the 
internal states. The occupation of the internal states can be read out by fluorescence 
imaging. 
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In Sec. 15.11 we first propose a measurement sequence for the cross-correlation 
signature. Then we discuss a proposal to directly measure the Hawking phonons in 
Sec. 15.21 The proposal in Sec. 15.11 allows to determine any part of the covariance 
matrix (see Eq. (JH])). Thus, it can be used to determine the emerging entanglement 
between the inside and the outside of the black hole (see Sec. I4.2.3H . Note that it has 
been proposed earlier how to detect entanglement in the motional degrees of ions [32J. 
The basic mechanism of all proposals is the coupling of the ion displacements to their 
internal levels with lasers. 



5.1. Measurement of Ion Displacements 

In this section we discuss how the cross-correlation signal analyzed in Sec. 14.2.11 can be 
detected in an experiment. We propose to measure correlations in the ion displacements 
by coupling the motional degrees of freedom of the ions to their internal states. 

First, we explain how to relate the momentum- momentum correlations we discussed 
(see Eq. (162]) ) to experimentally accessible correlations in the ion displacement and 
analyze how accurately the latter should be measured. To this aim we rewrite Eq. (1551) 
for a continuum system in a region of constant flow as 

8p(0) = ^-mc{e)n{6)^^-. (70) 

We can thus measure the momentum-momentum correlations by spatial derivatives of 
the ion displacements 

The analysis shown in Fig. H~2T b) confirms that we can use Eq. (1561) to get the order 
of magnitude of the momentum-momentum correlations for a finite ion ring. Thus, we 
can estimate the magnitude of the cross-correlation signal as 



'i+A 



^ j+A))H |yg^ <*y (72) 

1 L J m N A (d — Vi){Cj — Vj) 



where we used rii « nj w N/(2ir) and q ~ Cj ~ (2tt)/T. Therefore, the angle-angle 
correlations must be detected in an experiment with the accuracy 

e .= A { SW>,) = (yJ (73) 

In the following we are describing a setup to detect the angle-angle correlations of the 
ions with sufficient accuracy. We propose to illuminate the ion ring at two positions. 
The lasers should be focussed on one ion inside the supersonic and one ion inside the 
subsonic region on the ring. At each position a laser beam couples two internal levels \g) 
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and |e) of the ions with the transition energy coj. The lasers should fullfill the resonance 
conditions 

u = uj + kv— (74) 

for their frequencies and wavenumbers, which takes into account the Doppler shift kv. 
The ion traverses the pulse beam in the time T/N which is much shorter than the 
time scale T = 27r/w ro t of phonons at small wavenumbers. Thus, we can neglect the 
ion motion for the further analysis. After going to the frame rotating with uq and 
applying the rotating wave approximation the dipolar coupling Hamiltonian [ID] of one 
illuminated ion becomes 

n dip = hn (a + e^ s§ + F.c.) , (75) 

where Q is the Rabi frequency of the laser transition. Note that the ions remain in the 
Lamb-Dicke limit 

h 

k < 2tt (76) 



mNu mt 

during the experiment, i.e. kL58 <C Air 2 . 

We propose to prepare the internal state of each ion in the superposition (|g) + 
|e))/ y/2 before the experiment. After the creation of the black hole and the illumination 
with the lasers the probability that the two measured ions are in the states \gg) or |ee) 

p (56 u 56^ = i + (J^J sin 2 (2Qt) y (c*W (77) 

is measured through a repetition of the experiment. We propose to use the Rabi 
frequency 2Qt = tt/2. If the measurement is repeated M times, the standard deviation 
of the average number of binomially distributed events P is 



Af^ ^'-q ' , (7f 

1 M 2-JM 



It should be smaller than the required accuracy of the signal size and thus 



->(t)W M ^) ! 



2AT 3 (a -Vi){cj -Vj) 
A% 4 (k B T R /H) 2 



2 



(79) 



measurements are necessary. 

We are now calculating M for the example of N = 1000 9 Be ions, discussed 
in the beginning of this Sec. [5j We propose to use the A = 313nm transition in 
9 Be [5H]. We have checked that the cross-correlations remain clearly visible for the 
resolution A = 50 (see Fig. [T5|) . In this case M > 1.1 ■ 10 5 measurements are required. 
This basic measurement proposal can certainly be improved by employing additional 
techniques, e.g., using interferences with additional lasers. But we do not propose further 
experimental setups here, since these should be adapted to specific implementations of 
our proposal. 
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Figure 15. Correlations ((60, - 59 t+5 Q)(89 J - S0 ]+5O )) ■ m/(HT) ■ L 2 /(2ir) 2 at time 
t = 0.6T with real space lab frame positions. These correlations are analougous to the 
momentum-momentum correlations shown in Fig. |9Jd. We consider N = 1000 ions, 
crv min T = 2tt • 0.25, e 2 /47re = ^p-^-, and t = 0.05T. Full Coulomb interactions 
are considered. 



The accuracy of the classical equilibrium positions of the ions must satisfy |A#°| 2 < 
y/e, such that it does not influence the cross-correlation signal. Inaccurate equilibrium 
positions can be caused by inaccurate external forces (see |Appendix A[ ) during the black 
hole formation. Assuming that the external force is F{(t) = (1 + ^y)Fi(t) (see Eq. ( 1A.3I) ) 
with 7<1, the deviation in the classical positions fullfills the equation 



L d 2 A6i 
— m — — — 

2tt dt 2 



lF.it). 



Integrating Eg. [HOI we can estimate 



< 



7 



\(TV n 



dem , 

dt 

t 2 27T 



+ 



T 2 2vr |i^ c | 

2 N 2 



2 L 4ire m L 2 



- v„ 



7T 



For t = 0.05T we get |A0 4 | < 2tt • 0.257, thus the accuracy 7 < y/e/{2-K ■ 0.25) 
is required for the external forces. 
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5.2. Measurement of Hawking Phonons 

The measurement scheme described in the previous section requires the accurate control 
of the ion acceleration during the creation of the black hole (see Eq. (IHTj) ). Here we 
propose a scheme to avoid this difficulty (see also [H]). If the number of ions in the 
supersonic region (1 — 2a)N is small, the displacement of the ions due to the creation of 
the black hole will be small compared to the equilibrium motion of the ions. Then the 
creation of the black hole is adiabatic. In this case the measurement of cross-correlations 
is not possible. Instead, we describe in this subsection how to directly measure the 
emitted Hawking phonons outside of the black hole. 
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The following setup is studied: After the black hole formation in the small time 
interval r, Hawking phonons are emitted at the black hole horizon. We propose to detect 
these phonons by coupling the oscillation of the ions to their internal state with a laser. 
The laser drives a transition between two internal states |a) and |6) of the ions, which are 
prepared in the state | a) . It illuminates N ions outside of the black hole and is following 
the motion of these ions for the time interval t m . Since t m determines the spectral 
width of the laser, it must be large compared to the inverse Hawking temperature 
t m 3> h/ik^Tn). The laser frequency u should fullfill the resonance condition 



L 

u = tuj + co po + kv—, (82) 
2n 



which takes into account the transition frequency of the ions uj, the Doppler shift kv, 
and the relevant phonon frequency u po . After going to the frame rotating with ujj and 
u po and applying the rotating wave approximation the dipolar coupling Hamiltonian 
[441] in the Lamb-Dicke limit (see Eq. (!76|) ) becomes 



n dip = J2 (^<r}a P e*+ k W' N + H.c^j + £ h (u p - u po ) o) p a p , 



(83) 



where the sums extend over the illuminated ions j and the relevant phonon modes p. k 
denotes the wavenumber of the laser and 



2h 

Op = -ink a I (84) 

TTlUr, 



are the effective Rabi frequencies with the bare Rabi frequency Q. We describe here 
the coupling of the phonon modes to a sideband transition. In the limit of a laser 
pulse with small spectral width iV Nt m /T we can introduce the spin £ = N/2 
operator S + = exp (— ikLj/N) cr+/2 and transform it into a bosonic field with 
the Holstein-Primakoff transformation 



£+ = &V 2 £-&t&« bWN (85) 

For (6+6) w (ajop) < k B T n /{hw mt ) < (see Eq. (jST])), we can use £+ ~ \fNtf . The 
expectation value (tfb) gives the number of ions in the state \b). So the coupling of the 
relevant phonon modes to the illuminated ions can be described by the Hamiltonian 

Hdip = hVN^2 {^pb ] CLp + H.c) + ^ K^p - u Pa )o) p a p . (86) 
p p 

The number of excited ions after their illumination is 







12 
1 1 




^l 2 



^ M = E sin2 * /E i^' 2tm ( 8? ) 
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with the assumption of large Rabi frequencies V Ntt po ^> co po . Thus, the number of 
excited ions in state \b) is proportional to the number of phonons around the mode u po . 
The standard deviation of the phonon number measurement is approximately 



Now, we are presenting an example that satisfies the requirements of this 
measurement proposal. We consider N = 10 5 singly charged 9 Be ions with ion 
spacing L/N = 2fim and average rotation frequency u TOt = 2tt x 1.2kHz (e 2 /47re ~ 
0.7/ (2N)mL 3 /T 2 ). The black hole region contains (1 - 2a)N = 0.007iV and the 
horizon region 27iiV = 0.004iV ions. The angular velocity in the subsonic region is 
^min = 2it x 0.99/T and in the supersonic region f max = 2tt x 2.4/T. The Hawking 
temperature in this system is k^T^jfi ~ 106/T w 2ir x 20kHz. In the small time 
interval r = 0.05T of the black hole creation, the ions are normally traversing the angle 
A8 m 2it x 0.05, which is large compared to the size of the black hole region. For the 
measurement we propose to illuminate N = 200 ions for the time tu = T/A. In this 
case M = 100 repetitions of the experiment are sufficient. 

6. Conclusion 

In summary, we have discussed in this paper the details of a recent proposal to observe 
the Hawking effect with ions rotating on a ring [lj. We have described how to create 
an analog black hole spacetime in this system (see Sec. [2]). The horizon emits Hawking 
radiation with a thermal spectrum (see Sec. H]). We have analyzed the emergence of 
correlations and entanglement between the inside and the outside of a black hole after 
its creation (see Sec. 14. 2p . These correlations are a signature for the pair creation 
mechanism of Hawking radiation. 

In this paper we have deepened our analysis of the emerging entanglement (see Sec. 
14.2.31) . We cannot observe the creation of entanglement at too high initial temperatures. 
The generated cross-correlations, instead, remain present at arbitrarily high initial 
temperatures. Thus, we find the transition from the quantum to the classical Hawking 
effect (spontaneous versus stimulated emission). 

Nevertheless, current technology allows to measure the quantum Hawking effect 
in an experiment. We have presented a detailed discussion of realistic measurement 
techniques. It is possible to directly measure the cross-correlation signal (see Sec. 15. ip 
or the emitted Hawking phonons (see Sec. 15. 2p . 

To conclude, we expect ring traps to be extremely useful for future quantum 
simulations. They offer great opportunities, especially for studying translationally 
invariant systems. 
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Appendix A. External Forces 

In this Appendix we present the detailed form of the external forces appearing in the 
general Hamiltonian ([I]) and the harmonic Hamiltonian Eq. (EJ) . These forces are chosen 
to enforce the imposed equilibrium motion of the ions Bf(t) in Eq. ([2]). 
The Coulomb force on the zth ion tangential to the ring is given by 



(0 X (*),..., M*)) 



4vre 



J2F c m)-9 3 (t)) 



with 



F C (A9) = 7r 2 sign(sin(A#/2))- 



cos(A#/2) 



(a.i; 



(A.2) 



' L 2 sin(A#/2) 2 ' 

where 6i(t) gives the ion position at time t. The classical equations of motion are now 
mL/(2iT)6i(t) = F?(t) + F e (8i(t)). We determine the local external force F e from this 
equation such that it guarantees the imposed equilibrium trajectories of the ions 

L 



F e (6,t) 



+ 
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2ttT 2 
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9'U {9~\o)) 
dg'v 



2ttT 
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2tt 

d 9v B 



dv 
d 2 a 
dv 2 
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1 2 



d 2 V 



dv 

N-l 



(A.3) 



This force is time- independent if the parameter t> min is time- independent. 

In harmonic approximation (see Eq. (|5|) the Coulomb force and the external force 
are encoded in the force matrix T = (/y) 



(ssf) J)j(n J: ' (n ' ^• /r(/) - 



(A.4) 
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The contribution from the Coulomb force is 



Aixen 



'r WW- *?(*)) 



i = j 



with 



7T 



1 + cos(Afl/2)' 



sin(A#/2) 3 

The contribution from the diagonal external force is 
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Appendix B. Velocity Profile 



The classical equilibrium positions of the ions, thus their velocity profile, are imposed 
by the function g as specified in Eq. (J2]). We make the choice 



T 



f3 — ah 



x — l+a 
72 



< x < o — 7x 
—71 < x — o < 71 

o~ + < x < 1 — a — 72 
-72 < x - (1 - 0") < 72 

1 — (J + 72<^< 1 



(B.l) 



with a = (i; max - w min ) /2 and /3 = (v mSiX + w min ) /2 and h(s) = 15 /8s - 5 /4s 3 + 3/8s 5 . 
We use g(0) = to determine the equilibrium positions of the ions. In the supersonic 
region crt> min T < 6 < 2n — av min T the constant angular ion velocity is t> max (see Eq. 
PJ) and in the complementary subsonic region its constant value is v (9) = v min . Thus, 
the black hole horizon is located close to 6*h = o-v min T. This stepwise definition has 
the advantage that the velocity profile is flat inside and outside of the horizon, that the 
width of both regions can be adjusted, and that the width of the horizon regions can 
be adjusted. This flexibility is useful for the detection of correlation patterns in Sec. 
14.21 Also this profile must be sufficiently continuous to generate a physically allowed 
equilibrium motion. 
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In a part of this paper we dynamically create a black hole metric from a flat 
metric. We choose to reduce the velocity in the subsonic region from 2ir/T at t = 0, 
corresponding to homogeneously spaced ions, down to v m[n at t r according to 



2vr 

nun + ( "7^7 



exp 



(B.2) 



We choose a Gaussian profile to guarantee ^-(t = 0) = 0. In this case, the velocity 
profile becomes (compare with Eq. (j3j)) 



v(0, t) = ^-He)) + dg^ {9)) dv^ {t) 



T 



dv n 



dt 



(B-3) 



Appendix C. System Dynamics and Equilibrium State 

The quasi-free quantum dynamics of the harmonic system (jSJ) are governed by the 
classical linear equations of motion for the first (see Eq. (J7|)) and second moments (see 
Eq flS}). The equations for the first moments can be written 



0_ 

Ft 



(C.l) 



with the matrix Q = (Gi 



m 2 ? 



V L ) ' m 



The dynamics for the second moments are governed by the equation 

^T(t) = g(t).T(t) + T(t).g(tf. 



(C.2) 



(C.3) 



We determine the thermal state with temperature T of homogeneously spaced ions at 
rest by a mode decomposition of the system. The Fourier transform O diagonalizes the 
system © 



A? 



N 



56, 



to 



i=i 

JV 



i=l 



fc=l 



2^ + T^M^ 



(C.4) 



(C.5) 



The mode frequencies are 



fc.= (£)'(O/0V 



(C.6) 
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V m . n T/{2n) v mjn T/(27t) 



Figure Dl. log {[i — 1) for N = 100 ions, and av m i n T = 2tt ■ 0.25. In this stability 
diagram values of log (fi — 1) close to zero represent stable systems, (a) Only nearest- 
neighbor interactions considered; (b) Full Coulomb interactions considered. 



According to the Bose-Einstein statistic each mode is on average occupied by 

(«*> = 



cxp 



ksTo 



phonons. In equilibrium the first moments vanish = due to the parity symmetry 
of the harmonic Hamiltonian (J5J). The covariance matrix at temperature T is given by 



»(!) i^ol^o 
N ( 1\ 

(Spidpj) = hm^2 Oik^k { (n k ) + - J O 
k=i ^ ' 



kj) 

(SOiSpj) = 0. (C.8) 

Appendix D. Stability Analysis 

In this Appendix we describe a stability analysis of our system. After a brief introduction 
to the method, we show results of the numerical simulations and put them into the 
context of previous works [S]. Finally, we show that the observed instabilities are not 
important for our proposal. 

The solution of explicit linear differential equations like Eq. fIC.ip can be written 



3 

with the monodromy matrix U(t) satisfying the initial condition U(0) = 1 and the 
dynamics 

d t V(t) = G(t)-V(t). (D.2) 
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■ unstable system 






• stable system 













5 10 15 20 25 30 35 40 45 50 



t/T 

Figure D2. Largest eigenvalue of LP(t)U(f) at time steps t=nT. This quantity is 
an upper bound of observables. We consider N = 1000 ions, av min T = 2ir ■ 0.25, and 
t = 0.05T. We compare the stable system e 2 /47re = °-^ 53 (green dots) with the 
unstable system e 2 /47re = °2at 46 (bhie squares). Full Coulomb interactions are 
considered. 

If the system is periodic G(t + T) = G(t), Floquet theorem states 

U(t) = X(t) • e m (D.3) 

with the periodic matrix X(t + T) = X(t) and the constant matrix R (see e.g., [12]). 
Therefore, the stability of the motion is determined by the eigenvalue of the matrix 
U(T) = exp (RT) with the largest magnitude. We denote this magnitude by the 
dynamics are unstable for /i > 1. The monodromy matrix U does not only describe the 
first moments (see Eq. ( ID. II) ). but also the evolution of the covariance matrix (see Eq. 

(Q)) 

r(t) = u(t)r(o)u T (t) = (u(t) <g> u(*)) r(o). (d.4) 

In our case the periodicity with period T is not the highest symmetry. The system is 
also invariant under combined translations in time and space (see Eq. ( )46l) ). Thus U(T) 
follows from V(T/N) 

V(T) = (T-V(T/N)f , (D.5) 

where T is the index translation matrix. 

The results of the numerical stability analysis are shown in Fig. IDl( a) for nearest- 
neighbor interactions and in Fig. IDl( b) for full Coulomb interactions as a function of 
(cT/(2tt)) 2 = 2iVe7(47reo) ■ T 2 /(mL 3 ) and v min T/(2n). For av min T = 2ir x 0.25, is 
the mean of the squares of the two sound velocities (see Eq. (Tl6|) ) on the ion ring. The 
analysis shows a distinction between three parameter regions. If the system is supersonic 
on the whole ring (at small c 2 ) the system is mostly stable with stripes of instabilities. 
If the system is subsonic on the whole ring (at large c 2 ) the system is always stable. In 
the interesting intermediate regime a subsonic region coexists with a supersonic region 
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on the ion ring, here the stability analysis is most complex. In this case the system 
is unstable apart from stripes of stability that bunch up in the central region which 
is most interesting for experiments. So we find that the appearance of instabilities is 
closely related to the presence of a black hole horizon. The structure of the instabilities 
is very similar to the result of the stability analysis in [5]. Black hole laser instabilities 
[4~3] may contribute to these behavior. They occur in the presence of two horizon when 
particles bounce between the horizons and enhance themselves. 

This stability analysis detects exponential instabilities. We show now that these 
exponential instabilities are not important for the proposed experiment. The increase of 
the experimental quantities is bounded by the maximal magnitude v of the eigenvalues 

of (U^(T)) n (U(T)) n for n — 1, 2, v is shown in Fig. ID2l for an exponentially stable 

and an exponentially unstable system. The exponential instability becomes dominant 
for t > 5T. Therefore, it does not have a consequence for the proposed experiment 
performed during t <T. 

References 

[1] B. Horstmann, B. Reznik, S. Fagnocchi, and J.I. Cirac, Phys. Rev. Lett 104, 250403 (2010). 

[2] S.W. Hawking, Nature 248, 30 (1974). 

[3] W.G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 

[4] C. Barcelo, S. Liberati, and M. Visser, Living Rev. Rel. 8, 12 (2005). 

[5] L.J. Garay, J.R. Anglin, J.I. Cirac, and P. Zoller, Phys. Rev. Lett. 85, 4643 (2000); L.J. Garay, 
JR. Anglin, J.I. Cirac, and P. Zoller, Phys. Rev. A 63, 023611 (2001); P.O. Fedichev and U.R. 
Fischer, Phys. Rev. Lett. 91, 240407 (2003). 

[6] R. Balbinot, A. Fabbri, S. Fagnocchi, A. Recati, and I. Carusotto, Phys. Rev. A 78, 021603(R) 
(2008); I. Carusotto, S. Fagnocchi, A. Recati, R. Balbinot, and A. Fabbri, New J. Phys. 10, 
103001 (2008); J. Macher and R. Parentani, Phys. Rev. A 80, 043601 (2009). 

[7] S. Giovanazzi, Phys. Rev. Lett. 94, 061302 (2005). 

[8] T.A. Jacobson and G.E. Volovik, Phys. Rev. D 58, 064021 (1998). 

[9] U. Leonhardt and P. Piwnicki, Phys. Rev. Lett. 84, 822 (2000). 
[10] B. Reznik, Phys. Rev. D 62, 044044 (2000). 

[11] W.G. Unruh and R. Schiitzhold, Phys. Rev. D 68, 024008 (2003). 
[12] R. Schiitzhold and W.G. Unruh, Phys. Rev. Lett 95, 031301 (2005); 

[13] O. Lahav, A. Itah, A. Blumkin, C. Gordon, R. Shahar, A. Zaats, and J. Steinhauer, Phys. Rev. 
Lett. 105, 240401 (2010). 

[14] T.G. Philbin, C. Kuklewicz, S. Robertson, S. Hill, Friedrich Konig, U. Leonhardt, Science 319, 
1367 (2008). 

[15] G. Rousseaux, C. Mathis, P. Maissa, T.G. Philbin, and U. Leonhardt, New J. Phys. 10, 053015 
(2008); S. Weinfurtner, EW. Tedford, M.C.J. Penrice, W.G. Unruh, and G.A. Lawrence, Phys. 
Rev. Lett. 106, 021302 (2011). 

[16] F. Belgiorno, S.L. Cacciatori, M. Clerici, V. Gorini, G. Ortenzi, L. Rizzi, E. Rubino, V.G. Sala, 
and D. Faccio, Phys. Rev. Lett. 105, 203901 (2010). 

[17] T. Jacobson, Phys. Rev. D 44, 1731 (1991). 

[18] S. Corley and T. Jacobson, Phys. Rev. D. 57, 6269 (1998); T. Jacobson and D. Mattingly, ibid 
61, 024017 (1999). 

[19] W.G. Unruh and R. Schiitzhold, Phys. Rev. D 71, 024028 (2005); S. Corley, ibid. 57, 6280 (1998); 
R. Brout, S. Massar, R. Parentani, and P. Spindel, ibid. 52, 4559 (1995); Y.Himemoto and T. 
Tanaka, ibid. 61, 064004 (2000); H. Saida and M. Sakagami, ibid. 61, 084023 (2000). 



Hawking Radiation on an Ion Ring in the Quantum Regime 



42 



[20] R. Schiitzhold and W.G.Unruh, Phys. Rev. D 78, 041504 (2008). 
[21] W.G. Unruh, Phys. Rev. D 51, 2827 (1995). 

[22] J. Drees and W.Z. Paul, Z. Phys. 180, 340 (1964); D.A. Church, J. Appl. Phys. 40, 3127 (1969); 

B.I. Deutch, F.M. Jacobsen, L.H. Andersen, P. Hvelplund, H. Knudsen, M.H. Holzscheiter, M. 

Charlton, and G. Laricchia, Phys. Script. T22, 248 (1988). 
[23] G. Birkl, S. Kassner, and H. Walther, Nature 357, 310 (1992); I. Waki, S. Kassner, G. Birkl, and 

H.I. Walther, Phys. Rev. Lett. 68, 2007 (1992); T. Schatz, U. Schramm, and D. Habs, Nature 

412, 717 (2001). 

[24] D. Porras, F. Marquardt, J. von Delft, and J.I. Cirac, Phys. Rev. A 78, 010101(R) (2008). 
[25] R. Balbinot, A. Fabbri, S. Fagnocchi, and R. Parentani, Riv. Nuovo Cim. 28, 3 (2005). 
[26] A. Fabbri, I. Carusotto, R. Balbinot, and A. Recati, Eur. Phys. J. D 56, 391 (2010). 
[27] Quantum fields in curved spaces, N.D. Birrel and P.C.W. Davies, Cambridge University Press 
(1982). 

[28] R. Balbinot, S. Fagnocchi, A. Fabbri, and G.P. Procopio, Phys. Rev. Lett. 94, 161302 (2005); R. 

Balbinot, S. Fagnocchi, and A. Fabbri, Phys. Rev. D 71, 064019 (2005). 
[29] As already pointed out in [5] using these as Kruskal coordinates, one is able to recover the correct 

behavior of the density-density correlator in 1 dimension for a homogeneous BEC scaling as the 

inverse squared of the distance between the two points. 
[30] J. Macher and R. Parentani, Phys. Rev. D 79, 124008 (2009). 
[31] R. Schiitzhold and W.G. Unruh, Phys. Rev. D 81, 124033 (2010). 
[32] A. Retzker, J.I. Cirac, B. Reznik, Phys. Rev. Lett. 94, 050504 (2005). 
[33] C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schumacher, Phys. Rev A 53, 4 (1996). 
[34] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 (2002). 
[35] A. Botero and B. Reznik, Phys. Rev. A 67, 052311 (2003). 

[36] M.M. Wolf, G. Gicdke, O. Kriiger, R.F. Werner, and J.I. Cirac, Phys. Rev. A 69, 052320 (2004). 
[37] J. Eschner, G. Morigi, F. Kaler, and R. Blatt, J. Opt. Soc. Am. B 20, 5 (2003). 
[38] D. Hanneke, J. P. Home, J.D. Jost, J.M. Amini, D. Leibfried and D.J. Wineland, Nature Physics 
6, 13 (2010). 

[39] C. Monroe, D.M. Meekhof, B.E. King, S.R. Jcfferts, W.M. Itano, and D.J. Wineland, Phys. Rev. 
Lett. 75, 4011 (1995). 

[40] D.J. Wineland, C. Monroe, W.M. Itano, D. Leibfried, B.E. King, and D.M. Meekhof, J. Res. Natl. 

Inst. Stand. Technol. 103, 259 (1998). 
[41] R. Schiitzhold, Phys. Rev. Lett. 97, 190405 (2006). 
[42] G. Floquet, Ann. Ecolc Norm. Sup. 12, 47 (1883). 
[43] S. Corley and T. Jacobson, Phys. Rev. D 59, 124011 (1999). 



